The Turbulent Shock Origin of Proto— Stellar Cores 



Paolo PadoanQ 

Harvard- Smithsonian Center for Astrophysics, Cambridge, MA 02138 

Mika JuvelaQ, 

Helsinki University Observatory, Tahtitorninmdki, P.O.Box 14, SF-00014 University of Helsinki, Finland 

Alyssa A. GoodmanP] 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 

and Ake Nordlundf] 

Astronomical Observatory and Theoretical Astrophysics Center, 
Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark 



ABSTRACT 

The fragmentation of molecular clouds (MC) into proto-stellar cores is a central 
aspect of the process of star formation. Because of the turbulent nature of super-sonic 
motions in MCs, it has been suggested that dense structures such as filaments and clumps 
are formed by shocks in a turbulent flow. In this work we present strong evidence in 
favor of the turbulent origin of the fragmentation of MCs. 

The most generic result of turbulent fragmentation is that dense post shock gas traces 
a gas component with a smaller velocity dispersion than lower density gas, since shocks 
correspond to regions of converging flows, where the kinetic energy of the turbulent 
motion is dissipated. 

Using synthetic maps of spectra of molecular transitions, computed from the results 
of numerical simulations of super-sonic turbulence, we show that the dependence of 
velocity dispersion on gas density generates an observable relation between the rms 
velocity centroid and the integrated intensity (column density), cr(V ) — I, which is 
indeed found in the observational data. The comparison between the theoretical model 
(maps of synthetic 13 CO spectra), with 13 CO maps from the Perseus, Rosette and Taurus 
MC complexes, shows excellent agreement in the cr(Vo) — / relation. 

The a(Vo)—I relation of different observational maps with the same total rms velocity 
are remarkably similar, which is a strong indication of their origin from a very general 
property of the fluid equations, such as the turbulent fragmentation process. 

Subject headings: turbulence - ISM: kinematics and dynamics - individual (Perseus, 
Rosette, Taurus); radio astronomy: interstellar: lines 
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1. Introduction 

The importance of turbulence in the process of star 
formation was recognized long ago (Von Weizsacker 
1951), and was discussed in a seminal paper by Lar- 
son (1981). Several successive works have tried to use 
the observational data to relate different properties 
of MCs with the physics of laboratory and numer- 
ical turbulence, such as power spectra of kinetic en- 
ergy Q ' probability distributions of velocity and veloc- 
ity differences [| intermittency (Falgarone & Phillips 
1990; Falgarone & Puget 1995; Falgarone, Pineau Des 
Forets, & Roueff 1995), and self-similarity Q 

During the last decade, numerical simulations of 
transonic turbulence (Passot & Pouquet 1987; Passot, 
Pouquet & Woodward 1988; Leorat, Passot & Pou- 
quet 1990; Lee, Lele & Moin 1991 ; Porter, Pouquet 
& Woodward 1992 ; Kimura & Tosa 1993; Porter, 
Woodward & Pouquet 1994; Vazquez-Semadeni 1994; 
Passot, Vazquez-Semadeni & Pouquet 1995) and highly 
super-sonic turbulence (Passot, Vazquez-Semadeni, 
Pouquet 1995; Vazquez-Semadeni, Passot & Pouquet 
1996; Padoan & Nordlund 1997, 1999; Stone, Os- 
triker & Gammie 1998; MacLow et al. 1998; Ostriker, 
Gammie & Stone 1999; Padoan, Zweibel & Nordlund 
2000; Klessen, Heitsch & MacLow 2000), on relatively 
high-resolution numerical grids, have become avail- 
able, and very detailed comparisons between observa- 
tional data with turbulence models of MCs have been 
performed (Padoan, Jones & Nordlund 1997; Padoan 
et al. 1998; Padoan & Nordlund 1999; Padoan et al. 
1999; Rosolowsky et al. 1999; Padoan, Rosolowsky & 
Goodman 2000). 

Recent numerical studies of super-sonic magneto- 
hydrodynamic (MHD) turbulence have brought new 
understanding of the physics of turbulence. The most 
important results are: 

• Super-sonic turbulence decays in approximately 
one dynamical time, independent of the mag- 
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netic field strength (Padoan & Nordlund 1997, 
1999; MacLow et al. 1998; Stone, Ostriker & 
Gammie 1998; MacLow 1999). 

• The probability distribution of gas density in 
isothermal turbulence is well approximated by 
a Log-Normal distribution, whose standard de- 
viation is a function of the rms Mach number 
of the flow (Vazquez-Semadeni 1994; Padoan 
1995; Padoan, Jones & Nordlund 1997; Scalo 
et al. 1998; Passot & Vazquez-Semadeni 1998; 
Nordlund & Padoan 1999; Ostriker, Gammie & 
Stone 1999). 

• Super-sonic isothermal turbulence generates a 
complex system of shocks which fragment the 
gas very efficiently into high density sheets, fil- 
aments, and cores (this is the general result of 
any numerical simulation of super-sonic turbu- 
lence) . 

• Super-Alfvenic turbulence provides a good de- 
scription of the dynamics of MCs, and an ex- 
planation for the origin of dense cores with 
magnetic field strength consistent with Zeeman 
splitting observations. (Padoan & Nordlund 
1997, 1999)0 

We call "turbulent fragmentation" the process of 
generation of high density structures by turbulent 
shocks. Since random super-sonic motions are ubiq- 
uitous in MCs, turbulent fragmentation cannot be 
avoided: it is a direct consequence of the observa- 
tional evidence. In some analytical studies it is tac- 
itly assumed that turbulent fragmentation can be ne- 
glected if the kinetic energy of random motions is in 
rough equipartition with the magnetic energy. This 
assumption is wrong, because motions along magnetic 
field lines are unavoidable, and so turbulent fragmen- 
tation occurs via super-sonic compressions along the 
magnetic field lines, as discussed by Gammie & Os- 
triker (1996) and Padoan & Nordlund (1997, 1999). 

Numerical simulations of turbulence have been 
used to discuss the turbulent origin of MC structures 
by Passot & Pouquet (1987). Vazquez-Semadeni, 



This particular result is supported almost exclusively by our 
work. A significant fraction of the astrophysical community 
still favors the traditional idea that a rather strong magnetic 
field supports MCs against their gravitational collapse. An 
example of numerical work that favors the traditional idea of 
magnetic support is Ostriker, Gammie &; Stone (1999). 
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Passot & Pouquct (1996) and Ballcsteros-Paredes, 
Hartmann & Vazquez-Semadeni (1999) have used 
two dimensional simulations to argue that MCs are 
formed by turbulence. Padoan & Nordlund (1997, 
1999) have shown that super-sonic and super- Alfvenic 
turbulence can explain the origin of magnetized cores 
in MCs, including the observed field strength-density 
(B - n) relation (Myers & Goodman 1988; Fiebig & 
Giisten 1989; Crutcher 1999). The idea that proto- 
stellar cores and stars are formed in turbulent shocks 
has been previously discussed by Elmegreen (1993). 
In that work, an analysis of the gravitational instabil- 
ity of the cores can be found. Here, we present new 
strong observational evidence in favor of the turbulent 
shock origin of proto-stellar cores. Such evidence is 
based on the fact that dense post shock gas traces a 
gas component with a smaller velocity dispersion than 
lower density gas, since it maps regions of converging 
flows, where the kinetic energy of the turbulent mo- 
tion is dissipated. 

In §2 and 3, numerical simulations and observa- 
tional data, used in this work, are briefly described. 
In §4 we compute the rms flow velocity as a func- 
tion of the gas density, in simulations of super-sonic 
turbulence, and show that it decreases for increasing 
values of the gas density. In §5 we show that such 
general property generates an observable relation be- 
tween the rms velocity centroid and the integrated in- 
tensity (roughly proportional to the surface density), 
for the J=l-0 13 CO transition, and in §6 the same 
relation is found in the observational data. Results 
are discussed in §7, and conclusions are summarized 
in §8. 

2. Numerical Models 

The numerical models used in this work are based 
on the results of numerical simulations of super- 
Alfvenic and highly super-sonic MHD turbulence, run 
on a 128 3 computational mesh, with periodic bound- 
ary conditions. As in our previous works, the initial 
density and magnetic fields are uniform; the initial 
velocity is random, generated in Fourier space with 
power only on the large scale. We also apply an exter- 
nal random force, to drive the turbulence at a roughly 
constant rms Mach number of the flow. This force is 
generated in Fourier space, with power only on small 
wave numbers (1 < k < 2), as the initial velocity. The 
isothermal equation of state is used. Descriptions of 
the numerical code used to solve the MHD equations 



can be found in Galsgaard & Nordlund (1996); Nord- 
lund, Stein & Galsgaard (1996); Nordlund & Gals- 
gaard (1997); Padoan & Nordlund (1999). 

In this work we neglect the effect of self-gravity, 
although that can be described with a different ver- 
sion of our code (Padoan et al. (2000)). Here we 
compare the relative velocity of regions of MC com- 
plexes as a function of their gas density or column 
density. Such regions are distributed across the full 
extension of the MC complexes, that is several pc. In 
our numerical models, driven continuously by an ex- 
ternal force on the large scale, self-gravity is respon- 
sible for the collapse of gravitationally bound cores 
formed by the turbulent flow, but does not affect sig- 
nificantly the large scale flow. Since we assume that 
large scale motions in MC complexes are due to tur- 
bulence, and not to a gravitational collapse, we can 
neglect self-gravity. Similarly, we have neglected the 
effect of ambipolar drift, since it is not relevant for 
motions on the scale of several pc, although that is 
computed in a different version of our code, using the 
strong coupling approximation (see Padoan, Zweibel 
& Nordlund 2000). 

In order to scale the models to physical units, we 
use the following empirical Larson type relations, as 
in our previous works: 

/ T \ 05 

M - W \ja) • (1) 

where M. is the rms sonic Mach number of the flow, 
and a temperature T = 10 K is assumed, and 

<») = 2 '° x 103 (l^) • ( 2 ) 

where the gas density n is expressed in cm~ 3 . The 
rms sonic Mach number is an input parameter of the 
numerical simulations, and can be used to re-scale 
them to physical units. When comparing theoretical 
models with observations, M is in fact the only pa- 
rameter that we need to match (or its two dimensional 
equivalent on the maps, a v -see §3). In the absence of 
self-gravity and magnetic field, statistical properties 
of turbulent flows (very large Reynolds number) with 
the same value of M should be universal, and inde- 
pendent of the average density (for isothermal flows 
without self-gravity). If the magnetic field is present, 
the rms Alfvenic Mach number of the flow is also an 
input parameter of the numerical simulations (it de- 
termines the magnetic field strength). The physical 
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unit of velocity in the code is the isothermal speed of 
sound, C s , and the physical unit of the magnetic field 
is C s (4tt(p))^ (cgs). 

In this work we use two numerical models with 
rms velocity 3.4 and 1.7 km/s, which corresponds to 
M. k 13.0 and 6.5 respectively. Using the Larson type 
relations (1) and (2) we get L « 42 and 10.5 pc and 
(n) ~ 48 and 190 cm -3 . The one dimensional rms 
velocity for the two models is a v = 2.0 and 1.0 km/s, 
which are also recovered from the analysis of the syn- 
thetic spectral maps computed with these models. 
The values of a v have been chosen for the appropriate 
comparison with the observational data presented in 
§ 3.[] The magnetic field strength is B w 5 /iG in both 
models. 

Maps of synthetic spectra of molecular transitions 
are computed, using a non-LTE Monte Carlo radia- 
tive transfer code (Juvela 1997, 1998), from the three 
dimensional density and velocity fields generated in 
the numerical MHD experiments. The method of 
computing synthetic spectra was presented in Padoan 
et al. (1998). For the purpose of this work we use 
only one molecular transition, namely J=l-0 13 CO . 
Uniform temperature, T — 10 K is assumed for these 
radiative transfer calculations, in agreement with the 
isothermal equation of state used in the MHD calcu- 
lations. We are presently computing thermal equilib- 
rium models of MCs that we will use in a future work 
to study the effect of realistic temperature variations 
in molecular spectra. 

Since the spectral noise resulting from uncertain- 
ties in our radiative transfer calculations is always 
much smaller than the typical observational noise, the 
comparison of synthetic spectra with observational 
data can be done only after adding noise to the syn- 
thetic spectra, and the effect of noise on the statis- 
tical properties of the spectra need to be quantified 
(see §4). 

3. Observational Data 

We choose to use the J=l-0 13 CO transition be- 
cause it samples the range of values of column den- 
sity we are here interested in, and also because several 



9 The Taurus MC complex has a v ~ 1.0 km/s, and L ss 40 pc, 
while the Larson type relation (1) would give L as 12.6 pc. 
For the purpose of this work we are interested in comparing 
the observations with numerical models with similar rms Mach 
number, and we do not try to match the physical extension of 
each MC complex. 



large maps of molecular clouds are available in this 
transition. We compare maps of J=l-0 13 CO syn- 
thetic spectra with some of the largest observational 
J=l-0 13 CO spectral maps in the literature from the 
following Galactic regions: the Perseus MC complex 
(Billawala, Bally & Sutherland 1997), the Taurus MC 
complex (Mizuno et al. 1995), the Rosette MC com- 
plex (Blitz & Stark 1986; Heyer et al. 2000). These 
MC complexes have an extension of approximately 
30-50 pc, and radial velocity dispersions in the range 
1-2.4 km/s. We define the total rms velocity of a map 
as the rms velocity weighted with the total spectrum 
of the map, T(v): 



CTv = 



where 



_ J2v vT ( v ) dv 



(3) 



(4) 



The Blitz & Stark map of the Rosette MC complex 
has a v = 2.4 km/s, but, limited to the region that 
matches the more recent Heyer et al. map, the value 
is a v = 2.0 km/s (for both maps). The full map 
of the Perseus MC complex also yields a value of 
o v = 2.0 km/s. Taurus is instead much less turbulent, 
despite its large spatial extent, with a v — 1.0 km/s. 

The angular resolution is inversely proportional to 
the diameter of the antenna: 4 m for the Taurus map, 
7 m for the Perseus and the Blitz & Stark Rosette 
maps, and 14 m for the Heyer et al. map of Rosette. 
Assuming a distance of 140 pc for Taurus, 300 pc for 
Perseus, and 1600 pc for Rosette, the spatial reso- 
lution of the maps is 0.1 pc for Taurus, 0.15 pc for 
Perseus, 0.84 pc for the Blitz & Stark map of Rosette, 
and 0.42 pc for the Heyer et al. map of Rosette. The 
spectral resolution is 0.1 km/s for Taurus, 0.273 km/s 
for Perseus, 0.68 km/s for the Blitz & Stark map of 
Rosette, and 0.06 km/s for the Heyer et al. map of 
the same cloud. 

Rms noise N and average spectrum quality (Q) 
(see Padoan, Rosolowsky & Goodman 2000) also vary 
from map to map. The spectrum quality, Q, is related 
to the signal-to-noise, S/N. It is defined as the ratio 
of the rms signal (over the whole spectrum or inside 
a velocity window) to the rms noise, N: 



VE, T{vfdv 
N 



(5) 



The usual definition of S/N is based on Gaussian fits 
of the spectra, which we prefer to avoid because the 
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J=l-0 13 C0 transition typically yields spectra with 
significant non-Gaussian shape and multiple compo- 
nents. The spectrum quality is a sort of signal-to- 
noise weighted over the whole spectrum. The relation 
between Q and S/N is discussed in Padoan, Good- 
man & Rosolowsky (2000). We call average spectrum 
quality (Q) the value of Q averaged over the whole 
map. Values of {Q), N, resolutions and a v are listed 
in Table 1 for all the observational maps. In the fol- 
lowing, when different maps are compared with each 
other, we make noise and velocity resolution equal in 
the different maps, by adding noise and reducing the 
velocity resolution where necessary. 

4. Velocity Dispersion Versus Gas Density in 
Super Sonic Turbulence 

Turbulent fragmentation generates a complex sys- 
tem of dense post shock sheets, filaments and cores, 
reminiscent of structures observed in molecular cloud 
(MC) maps. An example of a two dimensional pro- 
jection of the three dimensional density field from a 
simulation of super-sonic turbulence is shown in Fig- 
ure 1. In previous works (Padoan, Jones & Nord- 
lund 1997; Padoan & Nordlund 1997, 1999), we have 
shown that, besides this morphological similarity, the 
density field of numerical super-sonic turbulence has 
important statistical properties in agreement with the 
density field of observed MCs. 

Figure 2 (left panel) is a two dimensional section 
(no projection) of the same three-dimensional density 
field used in Figure 1 . A complex system of filaments 
is apparent, with a number of high density cores in- 
side the filaments. In three dimensional super-sonic 
turbulence, filaments are formed by two dimensional 
compressions (at the intersection of sheets) and the 
densest cores are formed by three dimensional com- 
pressions. Most of the "filaments" in two dimensional 
sections like Figure 2, are instead two dimensional 
cuts through sheets, and most of the cores are local 
density maxima, due to fluctuations in the shock ve- 
locity (they usually corresponds to strongly curved 
segments of filaments) . Such density maxima are of- 
ten unstable to gravitational collapse, and are the ori- 
gin of proto-stcllar cores. 

Since the dense gas originates in shocks, that is in 
regions of converging flows, it should move with sig- 
nificantly lower velocity than the lower density tur- 
bulent flow. This is illustrated in the right panel of 
Figure 2, which shows the magnitude of the flow ve- 



locity on the same plane as the left panel. Dark blue 
is low velocity, and dark red high velocity. It is clear 
that high density filaments trace regions of low veloc- 
ity, at the intersections of high velocity "blobs" . This 
general property of super-sonic turbulence is quanti- 
fied by the dependence of the rms flow velocity, <r(v), 
on the gas density: 



a(v) = {(v - vf\p), 



(6) 



We repeat for different values of p, and than plot the 
result versus the gas density. Expression (6) means 
that the rms flow velocity is obtained as an average 
over the whole computational box, using only posi- 
tions where the gas density has a value contained in an 
interval centered around p. The average is repeated 
for different intervals of values of p, to span the whole 
range of densities. 

The plot is shown in Figure 3, at time t/tdyn = 0.07 
(square symbols) and t/td yn = 1.0 (asterisks), where 
t/tdyn is the time in units of the dynamical time (de- 
fined as the size of the computational box divided 
by the rms flow velocity), and the density is uni- 
form in the initial conditions, at t/td yn = 0.0. At 
t/tdyn = 1-0, regions of p « (p) have <j(v) comparable 
to the total rms velocity, a v — 3.0 km/s, while re- 
gions with 10 times higher density have much smaller 
rms velocity, a(v) w 2.3 km/sec, and at p = 100 (p) 
a(v) « 0.9 km/s. The rms velocity conditioned to 
the gas density decreases sharply with increasing gas 
density already at a very early time, t/tdyn — 0.07, 
when the density field is still very smooth. 

5. Rms Velocity Centroid versus Integrated 
Intensity 

Observational spectral maps of MCs, do not pro- 
vide a direct estimate of the three dimensional ve- 
locity field and gas density. Only radial velocity and 
velocity-integrated intensity (which is roughly pro- 
portional to the surface density) , averaged along each 
line of sight, arc directly available from the data. 
However, the dependence of the rms velocity cen- 
troid on the integrated intensity should resemble the 
a(v) — p relation, since lines of sight of high intensity 
are usually dominated by one or more dense cores. 
The velocity centroid, Vb, is the average velocity along 
an individual line of sight x on the map: 



V (x) 



J2 v vT(v,x)dv 



(7) 
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where T(v,x) is the signal (antenna temperature), at 
the velocity v and position x, and dv is the width of 
the velocity channel. The integrated intensity, J(x), 
is: 

7(x) = ^2>,x)<fo. (8) 

V 

We compute the rms velocity centroid, cr(Vb), con- 
ditioned on I: 

v(V ) = ((Vo-Vo) 2 \I) x , (9) 

(spatial average) and plot it against I. Expression 
(9) means that the rms velocity centroid is obtained 
as an average over the whole map, selecting all lines- 
of-sight on the map with values of I inside a given 
interval, [1,1 + dl]. This rms velocity centroid is 
not equivalent to a local line width, because it is ob- 
tained as an average over the whole map. The plot 
is shown in Figure 4, for two maps of J=l-0 13 CO 
synthetic spectra with different total line widths a v 
(defined in §3). Figure 4 shows that <j(Vq) decreases 
significantly for increasing values of I. The general 
property of super-sonic turbulence, namely high den- 
sity gas moves relatively slowly, is therefore apparent 
also in the observable relation <r(Vo) — I. The rela- 
tion c(Vo) — I is affected by noise and by the width 
of the velocity channels, which must be taken into 
account when comparing different maps. Noise has 
been added to the synthetic spectra used for com- 
puting the plot in Figure 4, to yield a value of the 
spectrum quality comparable to a typical value found 
in the observational data used in this work, (Q) = 3.5 
(see §3 for the definition of (Q)). 

In Figure 5 we show the effect of noise (left panel) 
and spectral resolution (right panel). The effect of 
increasing the noise (decreasing the value of (Q)) is 
that of making the cr(Vb) — / relation steeper, because 
the uncertainty in the determination of the velocity 
centroids due to noise contributes to the dispersion of 
velocity centroid values, and the effect increases at de- 
creasing values of integrated intensity I. Lower spec- 
tral resolution further increases the same effect. How- 
ever, we have verified that if the noise is low enough, 
(Q) > 20, both noise and velocity resolution have no 
effect on the <t(V ) — I relation, and the squared sym- 
bols in the left panel of Figure 5 correspond therefore 
to the intrinsic relations. We can conclude that this 
observable relation truly originates from the three di- 
mensional a(v) — p relation, with only a partial con- 
tribution from noise. 



We have also verified that spatial resolution does 
not affect our results. The spatial resolution can be 
decreased significantly, by rebinning the map to a 
smaller number of spectra, without any appreciable 
variation in the er(Vb) — I relation. As the spatial res- 
olution is decreased, however, the statistical sample 
(number of spectra) decreases, and statistical fluctua- 
tions (deviations around the high resolution a(Vo) — I 
relation) become progressively more important. 

6. Observational cr(Vo) I Relation 

We have shown in the previous section that the 
&(Vo) — I relation is sensitive to the value of the rms 
noise and to the spectral resolution. In the following 
plots, where we compare different spectral maps from 
observations and models, we have therefore added 
noise to the spectra and increased the velocity chan- 
nel width dv to match the map with the largest noise 
and dv. We have not modified the spatial resolution 
in any map, since that has no systematic effect on the 
cr(Vo) — / relation, as commented above. 

The left panel of Figure 6 shows the cr(Vo) — I re- 
lation for the maps of MC complexes introduced in 
§3. For the Rosette MC complex we have used the 
portion of the full Blitz & Stark map that matches 
the Heyer et al. map. The two models used for the 
comparison have <r v = 2.0 km/s, similar to Rosette 
and Perseus, and a v = 1.0 km/s, similar to Taurus. 
It is remarkable that the Rosette and the Perseus MC 
complexes have indistinguishable cr(Vo) — I relations, 
which are also coincident with the theoretical predic- 
tion (square symbols), for the same value of a v . The 
result for Taurus is also in good agreement with the 
a v = 1.0 km/s model. Horizontal shifts could be ex- 
pected in the plot, since different MC complexes can 
have different surface density. However, MCs and MC 
complexes are known to approximately follow the Lar- 
son relation between density and size (equation (2)), 
which implies roughly constant surface density, and 
small horizontal shifts in the plot. 

To compute the plots in the left panel of Figure 6, 
all maps have been treated to make their rms noise 
and velocity resolution equal. This is achieved by re- 
binning spectral profiles into a smaller number of ve- 
locity channels, and by adding noise, when necessary. 
In order to check that the cr(Vo) — / relations from 
different maps treated in this way are really compa- 
rable, we have computed the cr(Vb) — I relation for the 
Rosette MC complex using both the Heyer et al. map, 
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and the portion of the Blitz & Stark map that matches 
the region covered by the Heyer et al. map. The re- 
sult is plotted in the right panel of Figure 6. The 
velocity resolution of the Heyer et al. map has been 
decreased from dv — 0.06 km/s to dv — 0.68 km/s, 
and noise has been added, to match exactly the veloc- 
ity resolution and rms noise in the Blitz & Stark map. 
As can be seen in the right panel of Figure 6, after 
this drastic treatment of the higher resolution map, 
the er(Vo) — / relations for the two maps are prac- 
tically indistinguishable from each other, which sup- 
port the validity of the comparison of different maps 
(left panel of Figure 6. In the right panel, the case 
of the full Blitz & Stark map is also plotted (square 
symbols). The full map has a higher total rms ve- 
locity, a v = 1.4 km/s, than the portion that matches 
the Heyer et al. map, and its er(Vo) — / relation is 
therefore steeper, as expected. 

7. Discussion and Conclusions 

The origin of proto-stellar cores is a fundamen- 
tal problem in our understanding of the process of 
star formation. Models that have been proposed to 
describe proto-stellar cores are based on i) static or 
quasi-static equilibrium (e.g. Curry & Mckee 2000; 
Jason & Pudritz 2000), ii) thermal instability (e.g. 
Yoshii & Sabano 1980; Gildcn 1984; Graziani & Black 
1987), iii) gravitational instability through ambipolar 
diffusion (e.g. Basu & Mouschovias 1994; Nakamura, 
Hanawa, & Nakano 1995; Indebetouw & Zweibel 
2000; Ciolek & Basu 2000), iv) non-linear Alfvcn 
waves (e.g. Carlberg & Pudritz 1990; Elmegreen 
1990, 1997, 1999), v) clump collisions (e.g. Gilden 
1984; Kimura & Tosa 1996), vi) super-sonic turbu- 
lence (e.g Elmegreen 1993; Klessen, Burkert and Bate 
1998; Klessen, Heitsch, & Mac Low 2000; Padoan et 
al. 2000). 

Many detailed comparisons between observational 
data and models, which support the idea of the turbu- 
lent origin of the structure and kinematics of molec- 
ular clouds, have been presented in our previous pa- 
pers (Padoan, Jones & Nordlund 1997; Padoan et al. 
1998; Padoan et al. 1999; Padoan & Nordlund 1999; 
Padoan, Goodman & Rosolowsky 2000). Models of 
numerical turbulence can be compared with the ob- 
servations by computing i) synthetic stellar extinction 
measurements, ii) synthetic spectral maps of molec- 
ular transitions, iii) synthetic Zeeman splitting mea- 
surements, iv) synthetic polarization maps. 



Maps of synthetic spectra contain a lot of infor- 
mation about the kinematic, the structure and the 
thermal properties of molecular clouds, and can be 
analyzed with different statistical tools. In this work 
we have presented a new statistical method to an- 
alyze spectral-line maps, which is very useful be- 
cause it probes directly a very general property of 
super-sonic turbulence, that is the fact that dense gas 
traces a gas component with a smaller velocity disper- 
sion than lower density gas. This property arises be- 
cause the gas density is enhanced in regions where the 
large scale turbulent flow converges (compressions) 
and the kinetic energy of the turbulent flow is dis- 
sipated by shocks. If local compressions are instead 
due to local instabilities (e.g. gravitational instabil- 
ity, or gravitational instability mediated by ambipo- 
lar diffusion), and the large scale motions are only the 
consequence of local instabilities (as it should be in a 
self-consistent picture), gas density increases with the 
flow velocity dispersion (see for example the model by 
Indebetouw & Zweibel 2000), contrary to the obser- 
vational evidence presented in this work. 

The main conclusion of this work is that every 
model for the origin of molecular cloud structure 
and proto-stellar cores should be tested against the 
cr(Vb) — / relation. Turbulent fragmentation provides 
a realistic scenario for the origin of proto-stellar cores, 
which satisfies this new observational constraint. 

We are grateful to Edith Falgarone and Phil My- 
ers for discussions that stimulated this work, and to 
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TABLE AND FIGURE CAPTIONS: 



Table 1: Cloud name; maximum spatial extension; total rms velocity; telescope beam size; velocity channel 
width; rms noise; average spectrum quality (signal-to-noise); bibliographic reference. 

Figure |l|: Two dimensional projection of the three dimensional density field from a simulation of isothermal 
super-sonic turbulence with rms Mach number ~ 10. 

Figure ^: Left panel: Two dimensional section (no projection) of the same density field used for Figure 1. 
Most filaments are sections of sheets, and most dense cores are density maxima inside curved segments of filaments, 
formed by fluctuations in the shock velocity. Right panel: Modulus of the three dimensional flow velocity on the 
same two dimensional plane as in the left panel. Dark blue is low velocity, and dark red high velocity. The dense 
filaments on the left panel are commonly found in regions of small flow velocity, at the intersections of patches of 
high velocity. 

Figure Rms flow velocity, conditioned to gas density, versus the gas density, computed from a simulation 
of super-sonic turbulence with rms Mach number ~ 10. Squared symbols are for an early time, just 7% of the 
dynamical time after an initial condition with uniform density; asterisks are for a time equal to a dynamical time. 
Density values are binned over 20 logarithmic intervals between the average and the highest density. Within each 
interval, the density grows by a factor of 1.04 at the early time, and 1.27 at one dynamical time. 

Figure ^ Rms velocity centroid, conditioned to integrated intensity, versus the integrated intensity (see text 
for details). Two maps of synthetic spectra from MHD turbulence simulations are used, with different values of 
the total rms radial velocity (intensity weighted), a v = 1.0 km/s (diamonds), and a v = 2.0 km/s (triangles). The 
model with larger rms velocity has also a slightly larger maximum intensity (surface density) , although both models 
are rescaled into physical units using the Larson relations (see §2). Integrated intensity values are binned over 12 
linear intervals between the lowest and the highest values in the map. Within each interval, the density increment 
is ~ 1 K km/s (for the a v = 1.0 km/s model) and ss 2 K km/s (for the a v = 2.0 km/s model). 

Figure ^ Effect of noise and spectral resolution. Left panel: Conditioned rms velocity centroid versus intensity. 
Triangle and diamond symbols are the same as in Figure 4, while squared symbols have higher spectrum quality (Q) 
(signal-to-noise). Different values of (Q) are obtained by adding different levels of noise to the synthetic spectra. 
(Q) = 3.5 is typical of the observational maps used in this work. The plot is not sensitive to the decreasing noise, 
for (Q) > 20. Right panel: Same as left panel, but the squared symbols are for lower spectral resolution (the 
velocity channels of the synthetic spectra is increased from dv = 0.2 km/s to dv = 0.6 km/s). Clearly, both the 
level of noise and the spectral resolution must be taken into account when comparing different maps. Integrated in- 
tensity values are binned over 12 linear intervals between the lowest and the highest values in the map, as in Figure 4. 

Figure ^: Observational cr(Vo) — I relation. Left panel: Conditioned rms velocity centroid versus intensity for 
different MC complexes: Rosette, Perseus, and Taurus. Square symbols are for maps of synthetic spectra with total 
rms velocity a v = 1.0 km/s and 2.0 km/s. Right panel: Comparison of the two maps of the Rosette MC complex 
(see text for details). Integrated intensity values are binned over 10 linear intervals between the lowest and the 
highest values in the map. 
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MC 


L [pc] 


a v [km/s] 


dx [pc] 


dv [km/s] 


N[K] 


(Q) 


reference 


Taurus 


38 


1.0 


0.10 


0.10 


0.24 


2.3 


Mizuno et al. (1995) 


Perseus 
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0.15 


0.27 


0.24 
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Billawala et al. (1997) 


Rosette 
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Blitz & Stark (1986) 


Rosette 
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2.2 


0.84 


0.68 
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4.0 


Blitz & Stark (matching region) 
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2.0 


0.42 
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Heyer et al. (2000) 
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